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Abstract 

We introduce two explicit examples of polynomials orthogonal on the unit circle. Moments and the 
reflection coefficients are expressed in terms of Jacobi elliptic functions. We find explicit expression for these 
polynomials in terms of a new type of elliptic hypergeometric function. We show that obtained polynomials 
are orthogonal on the unit circle with respect to a dense point meausure, i.e. the spectrum consists from 
infinite number points of increase which are dense on the unit circle. We construct also corresponding explicit 
systems of polynomials orthogonal on the interval of the real axis with respect to a dense point measure. 
They can be considered as an elliptic generalization of the Askey- Wilson polynomials of a special type. 
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1. Introduction 



Let C be some linear functional defined on all possible monomials z n by the moments 

c n =£{z n }, n = 0,±l,±2... 



(1.1) 



(in general the moments c„ are arbitrary complex numbers). Then the functional C is defined on the space 
of generic Laurent polynomials V(z) = ^2^L_ Nl a n z " where a n are arbitrary complex numbers and N i 2 are 
arbitrary integers. Namely we have (by linearity of the functional) 

n=—N-i 

The monic Laurent biorthogonal polynomials P n {z) are defined by the determinantal formula [13] 



P n (z) = (A„)- x 



where A„ arc defined as 



c Ci 
C-l c 

Cl-n C2-r 
1 Z 



c c\ ... c„_i 
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Cl-n C 2 -n • • • C 

It is obvious from the definition (1.2) that the polynomials P n (z) satisfy the orthogonality property 

£{P n (z) Z - k } = h n S kn , < k < n, 
where the normalization constants are 

ho = Co, h n — A„ + i/A„. 
This orthogonality property can be rewritten in terms of biorthogonal relation [19], [13], 

£{P n (z)Q m (l/z)} = h n S nm 
where h n = A n+ i/A„ and the polynomials Q n {z) are defined by the formula 



Q n (z) = (A„)- x 



Co c-i 
Cl c 

Cn-1 C„_2 
1 Z 



C— n 
Cl-n 

C-l 

z n 



(1.2) 



(1.3) 



(1.4) 



(1.5) 



(1.6) 
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i.e. the polynomials Q n {z) are again Laurent biorthogonal polynomials with moments c„ = c_ n . 
Laurent biorthogonal polynomials satisfy three-term recurrence relation [13] 

P n +l(z) H~ d n P n (z) = z(P n (z) + (z)) (1.7) 

with some recurrence coefficients b n , d n - In fact, recurrence relation (1.7) uniquely determine Laurent biorthog- 
onal polynomials P n (z) under the standard initial conditions P_i =0, P a = 1. 

An important special case of the Laurent biorthogonal polynomials is obtained if 

c_„ = c„ (1.8) 

where c„ means complex conjugated moments with respect to c„. 

In this case the biorthogonal partners are Q n {z) — P n (z). If additionally we demand that 

A„ >0, n = 0,l,2,... 

then there exists a nondecreasing measure [i on the unit circle \z\ = 1 such that 

P n (e w )P m (e- w )d^(6) = h n S nm (1.9) 

Jo 

Corresponding polynomials arc called the Szego polynomials orthogonal on the unit circle. In what follows we 
will denote the Szego polynomials by <$> n (z) instead of P n (z). 
The Szego polynomials satisfy the recurrence relation 

* n+ i(*) = z* n (z) - a n z n $ n (l/z), (1.10) 

where a n = — $„ + i(0) are so-called reflection parameters (sometimes called Schur, Geronimus, Verblunsky. . . 
parameters). The normalization constants h n in the orthogonality relation (1.9) are expressed in terms of the 
reflection parameters as 

^ = (l-|«o| 2 )(l-K| 2 )...(l-|a„-i| 2 ) (l.H) 

If \a n \ ^ 1 then the orthogonal polynomials are nondegenerate, i.e. h n ^ 0, n = 0, 1, If however |ajv| = 1 for 

some N > 1 then we have a finite system of orthogonal polynomials on the unit circle: &o(z), &i(z), ■ ■ ■ &n(z). 
From the Szego recurrence relation (1.10) one can derive three term recurrence relation 

*„+i(«) + d n <P n (z) = z{$ n (z) + 6„$ n _i(«)), (1.12) 

where 

d n = b n = d n {\ - |a„_i| 2 ) 

a n -i 

Note that relation (1.12) is a special case of three-term recurrence relation (1.7) for the Laurent biorthogonal 
polynomials. 
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In what follows we restrict ourselves with the so-called real case, i.e. we will assume that c = 1 (the 
standard normalization condition) and all the moments c„ are real and symmetric c_„ = c„, whereas the 
reflection parameters a n are real and satisfy the restriction 

-Ka n <l, n = 0,l,... (1.13) 

In this case the polynomials < &„(z) have the orthogonality property 

,-2-k 

/ 1>n(e te )<3> m (e- te W(6) =h n S nm , (1.14) 
Jo 

where the measure possesses symmetric property 

p(2tt- 6) = -n{0) 
In case if the measure admits existence of the weight function p{6) 

dn(0) = p(9)d9 

we have that p{6) is a real nonnegative function on the unit circle (i.e. p{9) > 0, < 6 < 2ir) which is 
symmetric with respect to the reflection 

P (2tt -6)= p{6) 

Note that in the case (1.13) all expansion coefficients of the polynomials P n {z) are real, i.e. polynomials P n {z) 
take real values for real z. 

Orthogonality condition (1.14) in this case is equivalent to the conditions 

/ <5> n {e ie )e- lme dp{6) = 0, m = 0, 1, . . . n - 1 (1.15) 
Jo 

The recurrence relation for this case becomes 

* n+ i(z) - z<P n (z) - a n z n $ n (l/z), (1.16) 

with the initial condition $o( z ) = 1- Note that if the sequence of the real reflection parameters a n , n = 0, 1, 2, . . . 
is given then the monic polynomials $ n (z) are uniquely determined through the recurrence relation (1.16). Vice 
versa, if the real moments c„ are given (c = 1) with the symmetry condition c_„ = c„ and with the positivity 
condition A„ > 0, n = 1,2, . .. then all the reflection coefficients a„ are determined uniquely and satisfy the 

restriction — 1 < a n < 1, n = 0, 1, 

In what follows we will use a simple 

Proposition 1 Assume that —1 < a n < 1 /or a/i n = 0, 2, . . . and corresponding polynomials $ n (z) have the 
expression 

n 

$„(*) =^W„ s z s 

s=0 
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with some expansion coefficients W ns (W nn — 1 due to monicity of polynomials $„(z)J. Assume that there 
exists a real sequence c„, n = 0, 1, 2, . . . such that cq = 1 and 

n 

Y^W ns ~c s =0, n= 1,2,3,... (1.17) 
TTien £/ie sequence c n coincides with the (unique) sequence of moments c n corresponding to the polynomials 

Proof. Under condition — 1 < a n < 1 we have that the polynomials <& n {z) have a nondecreasing measure on 
the unit circle with some unique moments c„ (co = 1). The orthogonality relation for these polynomials can be 
presented in the form 

(£$„(z),z- m ) =0, m = 0,l,2,...,n-l (1.18) 

or in an equivalent form 

n 

W ns c s - m = 0, n= 1,2, m = 0, l,...,n- 1, (1-19) 

where c„ = (£, z n ) are the unique moments corresponding to the linear functional L. In particular, for m = 
we have the relation 

n 

Y,W ns c s = 0, n=l,2,..., (1.20) 

Now consider relation (1.17) for n = 1. We have c\ — aoc = 0. Hence c\ = a = ci because So = 1. We see that 
the moments c\ and c\ coincide. In order to prove that c„ = c„ for all n = 1, 2, 3, ... it is sufficient to apply 
induction with respect to condition (1.17). This finishes proving of the proposition. The importance of this 
proposition is that if we know the only orthogonality condition (1.17) then the general orthogonality relations 
(1.19) and (1.18) holds automatically and moreover the sequence c„ coincides with the moment sequence c„. 

Together with given polynomials <J> n (z) one can introduce the "reflected sign" polynomials $„(z) = (— 1)™$„(— z). 

Proposition 2 7/$„(z) are polynomials orthogonal on the unit circle with real reflection coefficients a n then 
the polynomials $„(z) are again polynomials orthogonal on the unit circle with the reflection coefficients a n = 
( — l) n+1 a n . The expansion coefficients for the polynomials $„(z) are W ns = (— l) n+s W ns . The moments c n 
corresponding to the polynomials l>„(z) are c n = (— l)™c„. 

The proof of this proposition is trivial. We will use this proposition in order to reduce the expansion coefficients 
W ns and reflection parameters a n to the most convenient form. 

Note also that for the case of real parameters a n we have simple expressions for $„(±1): 

n-l 

®n(l) = - a s ) (1.21) 

s=0 

and 

n-l 

$„(-!) = (-!)« JJ(l + (-l)«a a ) (1.22) 

s=0 

These formulas follow directly from recurrence relations (1.10). 
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2. Elliptic binomial coefficients and their properties 

Define the so-called "elliptic binomial coefficients" (EBC) E™. (In what follows we will often omit dependence 
of the Jacobi elliptic functions like sn(z; k) on the modulus k). If < j < n then 

J J. J. S n .I) s + 1 



and — otherwise. In (2.1) w is an arbitrary real parameter satisfying the restriction 

w ^ 4KM 1 /M 2 (2.2) 

with some nonzero integers Mi, M 2 . Recall that K = K(k) is the complete elliptic integral of the first kind. 

The elliptic modulus k is assumed to belong to the standard interval < k < 1. In this case all EBC are 
real for all j, n = 0, 1, 

We have the obvious symmetry property 

K-j=Ej, j = 0,l,...,n (2.3) 
Moreover, in the limit w = 0we obtain the ordinary binomial coefficients 

w^o 1 \jj j\{n-j)\ 
In the limit k = we get trigonometric degeneration of the elliptic Jacobi functions: 

fe^o J 1J- sm(w(s + 1)) 



2 Jl 



s=0 

where 



q = e~ Mw , (2.4) 



(2.5) 



are so-called the Gauss polynomials or q-binomial coefficients [3], [16] defined as 

nl _ (l- g ")(l-g"- 1 )...(l-g"--''+ 1 ) _ (q) n 
J\ q (l-q)(l-q*)...{l-ql) (</),•(</)»-/ 

where (a)„ = (1 — a)(l — aq) ... (1 — aq n ^ 1 ) is standard Pochhammer q-symbol (q-shifted factorial). 
Another limit k = 1 leads to hyperbolic functions degeneration. In this limit we have 

Hm EJ = - M» tMlp^i = (-1 ^7^^ , (2.6) 
fe-i 3 (q)j(q)n-j (-q)n (q)j(-q- n )j 

where q = e~ 2w . 

A less trivial property of the EBC is the recurrence relation: 

cn(wj)dn(w{n - j))E% = cn(wn)Ej^ + dn^n)^ 1 (2.7) 

In order to prove the relation (2.7) we need addition formulas for the Jacobi elliptic functions [24] from which 
the identity follows 

sn(wn)dn(w(n — j))cn(wj) = cn(wn)sn(wj) + dn(wn)sn(w(n — j)) 
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which is valid for all (complex) values of w, n, j. Using this identity we immediately arrive at recurrence relation 
(2.7). Note that in the limit w — relation (2.7) becomes the well known recurrence relation for the ordinary 
binomial coefficients: 

(M->cr) 

There are other recurrence relations for EBC of the similar type (2.7). For example if one replaces j — > n — j 
in (2.7) and apply the symmetry property (2.3) then one obtains 

cn(w(n - j))dn(wj)E? = cn^n)^" 1 + dn(wn)E^ (2.9) 

which differs from (2.7) by interchanging of the functions cn and dn. We present also the relation 

cn(w(n - j))E% = dn(w(n - j))E^ + cn^dn^j)^™- 1 (2.10) 

and the similar one 

cn(wj)EJ = dniw^E^- 1 + cn(wn)dn(w(n - j))E^ (2.11) 

obtained from (2.10) by the change j — > n — j. 
Another pair of relations 

dn(w(n - j))E] = cn(w(n - j))Ef^ + dn(wn)cn(wj) J B J "" 1 (2.12) 

and 

dn(wj)E r f = cniw^E^- 1 + dn(wn)cn(w(n - j))E r ^t (2.13) 

is formally obtained from (2.10) and (2.11) by interchanging cn and dn functions. 

All these relations are derived in a similar way using addition formulas for the elliptic Jacobi functions. 
Moreover, all these relations are reduced to the standard one (2.8) the limit w — > 0. 

Note that our definition of the "elliptic binomial coefficients" differs from the one given, e.g. in [11]. 

3. cn-polynomials on the unit circle 

Introduce the polynomials 

n 

& n C) (z) = J2WnsZ S , (3.1) 
s=0 

where the expansion coefficients W ns are defined as follows. If n = 0, 2, 4, ... is even then 

W na = (-I)* dn(w(n - s)) E^ (3.2) 

If n = 1, 3, 5, ... is odd then 

W ns = {-l) s+1 cn{w(n - s)) E^ (3.3) 

It is seen that W nn — 1 for all n — 1, 2, 3, . . . and hence the polynomials P n (z) = z n + 0(z n ~ 1 ) are monic. 
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Define the sequence a n = — $„ +1 (0), n = 0, 1, From definition (3.1) it follows that 

a 2n = cn(w(2n + 1)), a 2n +i = -dn(w(2n + 2)), ra = 0,1,2,... (3.4) 

For the reciprocal polynomials 3>Jj(z) = z n $ n (l/z) we have (by symmetry property of the elliptic binomial 
coefficients) similar expression 

n 
s=0 

for even values of n and 

n 
s=0 

for odd values of n. 

Proposition 3 The polynomials $1°^ (z) satisfy the recurrence relation 

$i%(z) = z<Pi c \z)-a n &%(z), n = 0,1,2,... (3.5) 

Indeed, comparison of terms z s leads to the conclusion that for the even values of n the recurrence relation 
(3.5) is equivalent to 

-cn(io(n + 1 - s))E r s l+1 + dn(w(n + 1 - s^E^ = -a n dn{ws)E™ 

where a n — cn(w(n+ 1)). But this relation coincides with (2.10). Quite analogously, for the odd values of n we 
see that condition of vanishing of all terms z s in expression (3.5) is equivalent to relation (2.12). 

From their definition it follows that the parameters a n satisfy an obvious restriction —1 < a n < 1 for all 
n = 0,1,.... Thus the polynomials Pi°\z) belong to the class of the polynomials orthogonal on the unit 
circle. From general properties [12] it follows that there exists a positive measure on the unit circle providing 
orthogonality property. 

We define a linear functional a by its moments 

c„ = (a, z n ) = cn(wn) (3.6) 

Proposition 4 If the functional a is defined as (3.6) then the polynomials ^^\z) satisfy the orthogonality 
property 

n 

(a,^ n c \z))=Y, W ^ = ^ « = 1.2,... (3-7) 

s=0 

Proof. If n is odd then the proof is trivial. Indeed, we have 

n n 

S = J2 W ns c s = - ]T(-l) s cn( W (n - s))cn(ws)E^. 

s=0 s=0 

Changing summation variable s — > n — s and using symmetry property (2.3) we have 

n 

S = -^2(-l) n - s cn(ws)cn(w(n - a))E% = -S, 

s=0 
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whence 5 = 0. 

Now consider the case when n is even. In this case we have 



S = J2 W ™ c ° = ^(-l) s dn(w;(n - s))cn(ws)E" s 



s=o s=a 

We can apply relation (2.7) to present the sum S in an equivalent form 

n n—1 

S = cn(wn) Y,(-iyE:-l + M^n) ^-l)^"" 1 

s=0 s=0 

But E™^ 1 = 0, hence we can replace s — > s + 1 in the first sum, then we have 

n-l 

S = (dn(ion) - cn(um)) ^(-l) 5 ^ 1 = 

s=0 

due to symmetry property (2.3) (note that in this case n— 1 is the odd and all terms in the sum vanish pairwise). 
Thus we proved the proposition for all values n = 1,2,.... 

Now apply Proposition (1) in order to conclude that condition (3.7) is equivalent to more general orthogo- 
nality condition 

(a,^ n c \z)z-i) =0, j = 0,1,..., n-l (3.8) 

or, equivalcntly to condition 

(a,<P^(z)<P^(l/z)) = h n 5 nm , (3.9) 

Moreover, we proved that the moments c„ providing orthogonality relation for the polynomials P^ c \z) are 
given by (3.6) and the functional a defined by (3.6) coincides with the orthogonality functional C. 

In order to get the orthogonality measure /z(0) on the unit circle corresponding to the moments c n = cn(wn) 
we use the Fourier expansion of the function cn(z) [24]: 

cn(z;fc) = -^ V exp^f* - l/2)/if ) (3 1Q) 

s= — oo 

where 

q = exp(-irK'/K) (3.11) 

and K' = K(k'), k' 2 = l-k 2 . 
We have 

Theorem 1 The polynomials <E>„(z) corresponding to the reflection parameters (3.4) satisfy the orthogonality 
property on the unit circle 

oo 

Ps$l?\z s )$£\l/z s ) = h n 5 nm , (3.12) 

s— — oo 

where 

z s = exp(mw(s-l/2)/K) (3.13) 
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and the discrete weights p s are given by 

IT 1 

P'=kK q s-l/2 + q s+l/2 ( 3 ' 14 ) 
Proof. It is sufficient to show that the measure given by (3.12) reproduces "true" moments c„, i.e. we should 
verify that 

oo 

c„ = cn(wn, k) = ^2 p s z™, n = 0, ±1, ±2, . . . (3.15) 

s— — oo 

But property (3.15) follows immediately from formula (3.10) if we put z = wn. Thus we verified that our 
measure gives correct values of all moments c„ = cn(wn). This finishes the proof, because in our case this 
measure is the unique due to the condition — 1 < a n < 1. 

In order to get expression for the normalization constants h n we use formulas (1.11) and (3.4) which yields 

n 

h n = M« Yl sn2 ( ws )i (3.16) 
s=l 

where [i n — k n for the even n and p n = fc™ _1 for the odd n. Clearly the normalization coefficients are all 
positive h n > 0, n = 0, 1, 2, . . .. 

Clearly for all s = 0, ±1,±2, . . . the spectral points (i.e. the only points of growth of the orthogonality 
measure) z s lie on the unit circle. Note that under the condition (2.2) we have that z s are dense on the unit 
circle. This means that for every point z s and for any small parameter e there exists another spectral point 
z s i such that \z s — z s >\ < e. In recent monograph by Barry Simon [20] one can find general results concerning 
existence of the singular continuous and dense point measures on the unit circle. Our cn-polynomials <J>„(z) 
provide a very simple explicit example leading to such measures. 

Sometimes the "reflected sign" polynomials & n C \z) = (— l) n ®n°\— z) arc more convenient. Using Propo- 
sition 2 we find that polynomials ^^\z) have the moments c„ = (-l)"cn(wn) and the reflection parameters 
a n = —cn(w(n + 1) for the even n and a n = —dn(w(n + 1)) for the odd n. The expansion coefficients are 
W ns — dn(w(n — s))E™ for even n and W ns = cn(w(n — s))E™ for the odd n. Transformation <&n°\z) —> &n\z) 
in this case is equivalent to the substitution w — > w + 2K. Indeed, we have c„ = cn(«m + 2Kn) = (— l)"cn(wn) 
due to the property cn{z + 2K) = — cn(z). This observation allows one to construct corresponding orthogonality 
relation for the polynomials & n C \z). We have almost the same relation as (3.12): 

oo 

]T Ps ¥ n c \z s )¥^\l/~z s ) = h n 5 nm , (3.17) 

s — — oo 

with the same weights p s and normalization coefficients and where 

z s = ex.p(inw(s — 1/2)/ K — in) — —z s (3.18) 



4. dn-elliptic polynomials 

In the previous section we considered the polynomials ^^\z) orthogonal on the unit circle with respect to the 
moments sequence c„ = cn(wn), n = 0, ±1, ±2, . . .. Consider now the polynomials ® n D \z) corresponding to the 
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moments given by the expression 

c„ = dn(wn), n = 0, ±1, ±2, . . . (4.1) 
Clearly, the moments (4.1) are real and satisfy the symmetry property c_„ = c n (it is assumed that < k < 

!)• 

As in the previous section it is verified that the reflection parameters a n corresponding the moments (4.1) 
are given by the expression 

a 2n = dn(w(2n + l);fc), a 2n +i = -cn(w(2n + 2); k), n = 0,1,2,... (4.2) 

The explicit expression of the polynomials $if\z) = X^"=o W ns z s in this case is very similar to expression 
for polynomials ^n\z). If n — 0, 2, 4, ... is even then 

W ns = {-l) s cn(w{n - s)) E™ (4.3) 

If n = 1, 3, 5, ... is odd then 

W ns = (-l) s+1 dn(w(n-s))EZ (4.4) 

Hence in this case the coefficients W ns are obtained from the coefficients (3.2) and (3.3) by simple interchanging 
the functions cn and dn. 

For real values of the parameter w the reflection parameters a n are real and satisfy the property —1 < a n < 1. 
If additionally the parameter w is chosen such that w ^ 2KM/N then we have — 1 < a n < 1 for all n = 0, 1, 2, . . . 
and hence the reflection parameters satisfy the condition (1.13) meaning that there exists a symmetric measure 
fi(8) providing orthogonality of the polynomials $if \z) on the unit circle. 

In a similar way as in the previous section we obtain the orthogonality measure 

Theorem 2 The polynomials <^^\z) corresponding to the reflection parameters (4-2) satisfy the orthogonality 
property on the unit circle 

oo 

Ps^ n D) (z s )^(l/z s ) = h n S nm , (4.5) 

s— — oo 

where 

z s — cxp(iTTws/K) (4.6) 

and the discrete weights p s are given by 

P* = 17 si -s ( 4J ) 

K q s + q s 

with q given by (3.11). 

The proof of this theorem follows immediately from the Fourier expansion of the Jacobi dn-function [24] 

dn(z;fc) = JL f expforz,/*) (48) 
v ' K ^ q s + q- s 



S — ~ OO 
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5. Elliptic derivative operator 

The polynomials <f> ( n\z) and $^ D) (z) are very close to one another in their explicit expression. There is a 
remarkable relation between them based on a notion of so-called "elliptic derivative operator" £. We define the 
operator £ on the space of all formal series YlsLo c s zS by its action on monomials: 

£z n = e n z n ~ 1 (5.1) 

where 

sn(wn) 
sn(uj) 

is the "elliptic number" . Note that in the limit ro = 0we have e n — n and the operator £ becomes the ordinary 
derivative operator: £ w= o = d z . In the limit k = we get the q-derivative operator 

n _ -n 
£z n = 1 1 z n-l 

q-q 

where q = exp(iw). 

The operator £ can be presented explicitly as an infinite sum of q-derivative operators. Indeed, introduce a 
set of q-derivative operators Wj which act by the formula 

WiM= f{ZeVe) ~/ z {Ze ~ Ve \ 3 = 1,2,3,... (5.2) 

with some fixed real parameter 9 (it is assumed that 9 ^ 2nM/N with integer M, N). On the set of monomials 
z n the operator Wj acts as 

Wjz n = sm(nj6)z n -\ n = 0,1,2,... (5.3) 

Hence the operator Wj sends any polynomial of degree n to a polynomial of degree n — 1 . 
Now we have 

oo 

£ = X>W 2j --i, (5-4) 

3=1 

where 9 — ttw/(2K) and 

Indeed, formula (5.4) follows from the Fourier expansion of the elliptic sn function [24]: 



, 2tt ^ sm((2j-l)t) 

^^lE gl /H_ r V 2 - t = *z/(2K) 

3=1 

The operator £ transforms any polynomial in z of degree n to a polynomial of degree n—1. 
We have 

Proposition 5 The polynomials $1°^ (z) and (z) are related as 



e n $i D _\(z)=£$W(z), e n ^} 1 {z)=£^\z) (5.5) 
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The proof is elementary and is based on explicit expression of the polynomials $£^(2), $> n D \z). 

As an elementary consequence of this proposition we obtain that the operator £ 2 transforms both families 
of polynomials $ n C \z) and $ n D \z) to themselves: 

£ 2 <P n c \z) = e n e n ^ n C 2 2 (z), £ 2 ^ n D \z) = e n e n .^ 2 {z) 



6. Limits k = and k = 1 

The "trigonometric" limit k = is not interesting because it leads to a degeneration. Indeed, in this limit, e.g. 
for the cn polynomials we have a2 n +i = — 1 f° r all n. Similarly, for dn polynomials in this limit we have a 2n = 1- 
This means that we have the degenerated case of the polynomials orthogonal on the unit circle. It is instructive 
to consider what happens with orthogonality relation in these degeneration cases. In the limit k = we have 
[24] K(0) = 7r/2 and K' = log(4/fc) + 0{k 2 log(fc)). Hence q -> k 2 /16 -> 0. For the case of cn-polynomials we 
have that only terms po and p\ in formula (3.14) remain nonzero in the limit k = 0. Indeed, in this limit we 
have only two spectral points z = e~ lw and z\ — e lw having two equal concentrated masses po = Pi = 1/2. For 
the case of dn-polynomials we have the only spectral point Zq = 1 with the concentrated mass p = 1. Thus in 
the trigonometric limit the problem becomes trivial. 

Consider the "hyperbolic" limit k = 1. In this case the reflection parameters a n for both cn and dn 
polynomials coincide: 

(-D™ 
cosh(w(n + 1)) 

The moments for both polynomials are 

c„ = 1/ cosh(wn) (6.2) 

Hence both polynomials coincide in this limit and we have from (2.6) the following expression for the power 
coefficients 

( _ ir+Sg -n/2 +s /2 {q -n )s[ _ q)s 

Wns = 2 r~; — ~ 7T7 — r^rr *- 6 - 3 ) 

1 + 5 ™ (3)4-9 n )s 

We can present polynomials <& n {z) in the form 

*»W - E = ( ? rl7; -*?" 1/a ) , (6-4) 

where 2 ^i stands for the basic hypergeometric function defined as [16] 

(clearly if a = q~ n with positive integer n then we have a finite sum of n + 1 terms). 

The reflection parameters a n satisfy the condition — 1 < a n < 1, n = 0, 1, . . ., hence there exists a unique 
weight function p{9) on the unit circle such that 

c„ = 1/ cosh(wn) = / p(0) cos(n0)d9, n = 0, 1, 2, . . . (6.5) 
Jo 
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Proposition 6 The weight function p(9) for the moments (6.2) is expressed in terms of the Jacobi dn function: 

p{6) = ^ dn(K(k)9/n), (6.6) 

where K(k) is complete elliptic integral of the first kind and the elliptic modulus k is determined from the 
transcendent equation 

w = irK'(k)/K(k) (6.7) 

Proof of this statement immediately follows from the Fourier expansion (4.8) of the dn function. The function 
p{9) satisfies the required symmetry property p(2tt — 9) = p{6) which follows from the periodicity property 
dn(u + 2K) = dn(w). Note that to every value of w > there corresponds a unique value of the modulus k in 
the canonical interval (0, 1) such that k — > 1 when w — > and k — > when w — > oo. 

The obtained polynomials P n {z) coincide with a special case of the family of biorthogonal polynomials 
introduced by Pastro [19]. Indeed, the Pastro polynomials depend on two parameters a, (3 and are expressed as 

P n (z)=A n2 ^^},zq^ 

where A n is a constant not depending on z. Comparing with our expression (6.4) we find that q"- 1 = q 13 ' 1 = -1 
(we need also to rescale the argument z — > —z/q). Thus our polynomials correspond to a special case of the 
Pastro polynomials with fixed values of the parameters a, (3. It is interesting to note that the weight function 
p{6) is expressed in terms of the elliptic function. This can be also derived from the results of Pastro [19] who 
obtained explicit expression of the weight function p{9) for arbitrary values of the parameters a, (3. Pastro 
presented his result for p(9) in terms of an infinite product. In our case this infinite product coincides with 
corresponding formula for the function dn(x). 

The "reflected sign" polynomials & n (z) = (— l)™<E> n (— z) have the moments c„ = (—1)"/ cosh(wn). It is 
easily seen that corresponding weight function on the unit circle for these polynomials is 

n k'Klk) 1 

7. Special case of polynomials orthogonal on regular polygons 

So far, we assumed that the parameter w satisfies the restriction w ^ AKMi/M 2 with some integer M 1 ,M 2 - 
This restriction is necessary in order for polynomials $> n (z) to be nondegenerate. 

In this section we consider the case when this restriction is omitted. This means that a degeneration occurs, 
i.e. ajv_i = ±1 for some positive integer N. Assume that under this condition the polynomial $jy(z) has only 
simple zeros: 

$at(z) = (z - z ){z - Zl) . . . (z - Z N -l) (7.1) 

i.e. Zi ^ z k when i ^ k. It can be easily showed [12] that under conditions \a n \ < 1, n = 0, 1, . . . , N— 2 all these 
zeros z s lie on the unit circle: \z.\ = 1. 
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In this case we can consider only a finite set of polynomials $o(z), &i(z),. . . &n-i(z) which are orthogonal 
on the finite set of these zeros on the unit circle [12] 

JV-l 

Ps^n{zs)^m{^/z s ) = K 6 nm , n, m = 0, 1, . . . , N - 1 (7.2) 

s=0 

where the discrete weights p s can be expressed as [26] 

_ h N-l (7 

Ps ~ ^(l/z s W N (z s ) (i - 6) 
As an example, consider the case of cn polynomials and condition 

w = KM/N, M = 1,3,..., N -2 (7.4) 

(it is assumed that M is co-prime with N). From (3.4) we see that under this condition a^N—i = — 1. Thus we 
need to calculate zeros of the polynomial &2n(z). We cannot substitute directly the value (7.4) into the explicit 
expression for $>2n(z) because in this case there is indeterminacy in expansion coefficients. Nevertheless we 
can put w — KM /N + e and then take the limit e = (because from recurrence relation (3.5) it is seen that 
expansion coefficients of the polynomials <& n {z) depend on w continuously). Then we easily obtain that under 
condition (7.4) only the last and the first terms survive in the power expansion of $>2n(z): 

$2N{Z) = Z 2N + 1 

and hence all zeros are simple: 

z s = exp(7rz(s- l/2)/N), s = -N + l,-N+l,...N (7.5) 

Note that these zeros coincide with corresponding spectral points (3.13) when to = KM/N and M is co-prime 
with N. 

We thus obtained explicitly spectral points z s in the finite-dimensional case. In order to get expression for 
the weights p s we should calculate the expression 

2JV-1 

*2JV-i(l/««) = (-l) m+1 cn(if (27V - m - l)M/N)E^~' i cxp(-i7rm(s - 1/2) /N) (7.6) 

m=0 

However we can avoid direct calculation of the sum (7.6) if we can construct a finite discrete measure on the 
unit circle which gives us appropriate values for the moments. In what follows we will consider only the simplest 
case when N is prime number and M = 1. We have the set of 2N moments 

c n = cn(Kn/N), n = 0, 1, . . . , 2N - 1 

We should find such values pj that 

JV N 

c n = cn(Kn/N)= ^ Pj z] = ^ Pj exp(Mrn(j - l/2)/N) (7.7) 

j = -N+l j=-N+l 
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To do this we return to the Fourier series (3.10) for the cn function and put z — Kn/N. It is convenient to 
present the summation variable s modulo 2N, i.e. s = j + 2Nm, where j = —N + 1, — N + 2,...N and m can 
takes all integer values — oo < m < oo . Then formula (3.10) can be presented in the form 



N N 



cn(Kn/N) = ^- ]T S(j; N) = JL ]T zj S(j;N) (7.8) 



kK ^ w ' ' kK 

j=-N+l j=-N+l 



where 

qj+2Nm-l/2 _|_ q-j-2Nm+l/2 



(7.9) 



m— — oo 

Wc thus have 



Pi = -£ R SV;N) (7.10) 



In order to calculate the sum S(j;N) explicitly we need the 
Lemma 1 For any q such that < q < 1 and any real a introduce the function 



oo 

F(a;q)= Y — . (7.11) 



S — — OQ 



Then the function F(a; q) has two equivalent representations: either 

1 (_ Q l+2a _ l-2a Q 2 2 2\ 

or 

F(a; g) = — dn(2aK'; k'), (7.13) 

7T 

where as usual q = cxp(—nK'/K). We adopt standard notation for the q-shifted factorial [11] (z;q) n = (1 — 
z)(l — zq) ... (1 — zq n ~ r ) and (z\, z 2l ■ ■ ■ , zjy; q) n stands for product of q-shifted factorials {z\; q) n . . . (zjy; q) n . 

The proof of this Lemma follows directly form the Ramanujan summation formula for the bilateral iV'i basic 
hypergeometric function [11]: 

^ {aq;q) n (b, aq, aq/c, bq/c; q) x 



s— — oo 



where (x; g) n means q-shifted factorial. In our case it is sufficient to put c = ab, a = —q a , b = q 1 / 2 . Then the 
Ramanujan formula (7.14) gives 

Ffcro 1/2 ) = V - = - [ —^- '—2- ' g ' g ' q) °° 

^ ^ ' s ^q(n+a)/2 + q -(n+a)/2 g a/2 + q - a /2 (_ g l+a _ g l-a >g l/2 ) q l/2. ^ 

Then we replace q —> q 2 and obtain formula (7.12). Formula (7.13) is then obtained from the well-known 
representation of the Jacobi elliptic functions in terms of infinite products [24]. Another method to prove (7.13) 
is using the well known Poisson summation formula [2] 

00 1 00 / 97rs\ 

f(t + sT) = - f[ — )oM^st/T) (7.15) 
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where the functions /(#), f(x) are connected by the Fourier transform: 

/oo 
f(t)e- M dt 
-oo 

Applying formula (7.15) to the function f{x) — 1/ cosh(/3a;) with q — e _/3 and using the Fourier expansion (4.£ 
of the dn- function we arrive at formula (7.13). 
Now we see from (7.9) and (7.12) that 

s(r,N) = F {^^; q 2N 

and we can express the sum S(j; N) in the one of two equivalent forms (7.12) or (7.13). In the "product form" 
(7.12) we have 

1 (_ n 2N+2j-l _ 2N-2J+1 AN AN . AN \ 

U ' '~ gJ-V2 + gV2-J (-g4JV+2j-l ; _ q iN-2j+l^ q 2N^ q 2N. ^JV)^ V - 10 ) 



S(j;A0 = - dn (i±lilk'-~k'\ (7.17) 



In the "dn-form" we have 

,vri- ATI = ... 

7T V N 

where K, K' and k, k correspond to q = q 2N = cxp(—2nNK'/K) instead of q. Such transformation (from q to 
q 2N ) correspond to so-called main 2iV-order transformation of Jacobi elliptic functions [1]. In more details, we 
have [1] 

where 



11 sn 2(lf;fc) 

It is interesting to note that the dn-function in (7.17) can be explicitly expressed as a rational function 

S(j;N) = R N (4>(j)) 

of the function 

We will not consider this form here (see, e.g. [1] for details, where relation with so-called Zolotarev rational 
functions of the best approximation is discussed). 
We thus have the following 

Proposition 7 Assume that w = K/N , where N = 1, 2, 3, . . . . Define the reflection parameters 

a-2n — cn(w(2n + 1)), a 2n+ i = -dn(w(2n + 2)), n = 0,1, . . . , N - 1 
such that a27v-i = — 1- Define the monic polynomials Q n (z),n = 0, 1, . . .2N by recurrence relation (1.16). 
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Then the polynomial &2n{z) — z 2N +l has simple zeros z s — cxp(ni(s—l/2)/N), s = —N+l, —N+l, . . . N. 
The polynomials $ n (z)> n — 0,1, ... , 2N — 1 are orthogonal on the finite set z s of spectral points on the regular 
2N-gon: 

N 

Pj®n{z j )§ m {l/zj) = h n 8 nm (7.19) 

j=-N+l 

where the discrete weights pj are given by formula (7.10) with S(j;N) given by (7.16) or (7.17). 

Remark. One can consider the generic case w — KM/N with arbitrary co-prime integers M,N. This will lead 
to slightly modified formulas for the weights p s . The case of dn-circle polynomials can be considered in a similar 
way. 



8. Polynomials orthogonal on the interval 

Now we can relate with our polynomials orthogonal on the unit circle P n (z) some polynomials S n (x) orthogonal 
on an interval of the real axis. This relation was first proposed by Delsarte and Genin [7] (see also [25] for 
further results concerning the Delsarte-Genin transformation (DGT)). 

Recall the main properties of the DGT. Given polynomials P n {z) with real reflection parameters a n satisfying 
the restriction — 1 < a n < 1 we construct the polynomials S n (x) by the formula 

o (T] _ Z- n ^ n (z) + Z n ^ n (l/z) 

Sn{x) 2»(l-a B _ 1 ) ' (8 ' 1} 

where x — (z 1 / 2 + z -1 / 2 )/2 and it is assumed that z ±1 / 2 = r 1 / 2 e ±* e / 2 when z = re l8 . It is easily verified that the 
polynomials S n (x) are monic n-th degree polynomials in x: S n (x) = x n + 0(a; Il_1 ). Moreover, the polynomials 
S n {x) are symmetric, i.e. S n {— x) = (— l) n S n (x) and satisfy the recurrence relation 

S n+1 (x) + VnSn-^x) = xS n (x) (8.2) 



where 



v n = j(l + a„_i)(l - a„_ 2 ) (8.3) 



Note that for n = we have formally w = due to assumed condition a_i = —1. Due to the same condition 
we have also v\ = (1 + a )/2. 

Under restriction — 1 < a n < 1 the recurrence coefficients v n are strictly positive u„>0,n=l,2,.... Thus 
the polynomials S n (x) are orthogonal on the interval [—1,1] with respect to a positive weight function 

J S n (x)S m (x)w(x)dx = K n 5 nm (8.4) 



where 



W { X ) = ^y * = 006(^/2) (8.5) 
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and 

Kn = VlV 2 ...V n = 2 1 - 2n h n {l - a^)" 1 

The orthogonality relation can be presented in an equivalent form 



J w{x) 



S n {x)S m (x) + S n (—x)S m ( x) 

dx = k„ d nm (8.6) 

which follows from the property w(— x) = w(x) for all symmetric polynomials S n (x). Relation (8.6) is more 
convenient when the measure is discrete one. 
Introduce the moments 

M n = j w(x)x n dx, n = 0,1,2,... 

corresponding to the weight function w(x) of the polynomials S n (x). It is easily seen that the moments M n are 
expressed in terms of the moments c„ as [7] 

M n = 2- n J2(%i-n/2 (8-7) 

j=0 

when n is even and M n = when n is odd. 

In the case of cn-polynomials we obtain that corresponding symmetric polynomials S n (x) are orthogonal 

s— — oo 

on the dense set of the points 



x s = cos 



TTW{S - 1/2) 



2K 

with the weights p s is given by (3.14). 

In the case of dn-polynomials we obtain similar orthogonality relation (8.8) but now the orthogonality grid 

is 

f1TWS\ 

and the discrete weights p s are given by (4.7). 

As expected, both cn and (in-type polynomials S n (x) are orthogonal with respect to measures which are 
positive and dense on the interval [—1,1]. 

It is well known [6] that to any family of symmetric orthogonal polynomials S n (x) one can associate two 
families of orthogonal polynomials P n (x) and Q n (x) by the formulas S2n(x) = 2~ n P n (2x 2 — 1) and S^n+i^) = 
2- n xQ n (2x 2 - 1) . Polynomials P n (y) and Q n (y) are monic P n (y) = y n + 0(y n - 1 ), Q n (y) = y n + 0(y n - 1 ) and 
orthogonal on the same interval [—1, 1]: 

«2n S nm (8.9) 

and 

J i Q n {y)Q m {y)wW{y + l)/2){y + lfl 2 dy = 2 2n ^ n 2n+1 S nm (8.10) 
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The polynomials P n {x) satisfy the recurrence relation [6] 

P n+1 {x)+b (x) = xP n (x), (8.11) 

where 

U n = 4u 2 nU2n-lj &n = 2(u 2 „ + V 2n +l) - 1 (8.12) 

Analogously, polynomials Q n (x) satisfy the recurrence relation 

Q n +i(x) + b n Q n (x) + u n Q n -\{x) = xQ n (x), (8.13) 

where 

U n = 4u 2n «2n+l, = 2(u 2ll+2 + "2n+l) - 1 (8-14) 

Using the representation (8.3) we obtain explicit formulas for the recurrence coefficients of the polynomials 
Pn(x): 

U n = -(1 + a 2 „-l)(l - aL- 2 )( 1 ^ a 2n-3), &n = ^ ( a 2n(l - 02n-l) - <l2n-2(l + 02n-l)) (8-15) 

Similarly, for the recurrence coefficients of the polynomials Q n (x) we have 

«n = |(1 + 02n)(l - oL-iX 1 - fl 2n- 2 ), &n = ^ (a 2 n+i(l - a 2n ) - a 2 „-i(l + a 2 „)) (8.16) 

Formulas (8.15), (8.16) were first obtained by Geronimus [12] (see also [20]). 
Explicitly polynomials P n (x) are expressed in terms of § n {z) as 

F -'"° '"*r.i;' + :" ti :" / '' *-e+o/> <»."> 

z"(l - a 2n _i J 

(formula (8.17) is due to Szego). Similarly for polynomials Q n {x) we have 

Qn(*)- 2 „ (zl/2 + z -l/ 2 )(l- a2n ) ' X={Z + Z ^ 2 ( 8 - 18 ) 

In what follows we will consider only the polynomials P n (x). All formulas for the companion polynomials 
Qn{x) are obtained in a similar way. 

Polynomials P n (x) are orthogonal on the interval [—1,1] with respect to the weight function 

V ; |sin0| VT^2 
Assume that polynomials <& n {z) have the expression 

n 

<P n (z) = J2w ns z s 

s=0 

with some coefficients W ns . Then from (8.17) we obtain corresponding expansion formula for polynomials P n (x) 
[5]: 

P n (x) = 2 1 -"(1 - a^-r)- 1 ^W 2n ,„ + J2(W 2n , n+s + W 2n , n - s )T s {x)^ , (8.20) 
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where T n (x) = cos(#n) is the Chebyshev polynomial of the first kind (it is assumed that x — cos 9). 

Thus if polynomials <& n {z) on the unit circle have explicit expression as expansion in terms of monomials z n 
then the corresponding polynomials P n {x) on the interval have explicit expression as expansion in terms of the 
Chebyshev polynomials T n {x). 

Return to our elliptic cn-polynomials. Corresponding polynomials P n {x) on the interval have the explicit 
recurrence coefficients 

sn 2 (>(2n - 1))(1 - dn(2tun)(l + dn(2w(n - 1)) 



u„ = 



; 

cn(w(2n + 1))(1 + dn(2wn)) - cn(w(2n - 1))(1 - dn(2wn)) 



2 

Using explicit expression (3.2) for the expansion coefficients W ng we obtain in case of the cn polynomials 



P n {x) = (-1)"2 1 -"(1 + dn(2wn))- 1 \ dn(wn)E* n + ^(-l) s ^« s ( d n(w;(n - a)) + dn(w(n + s))T s (x)j 

(8.21) 

These polynomials satisfy the orthogonality relation 

oo 

^ p s P»(x s )P m (x s ) = H n S nm , H n = uiu 2 ...u n (8.22) 

s — — oo 

where the spectral points are 



x s = cos 



7TW(S - 1/2) 



K 

and the weights p s are given by (3.14). 

The polynomials Q n {x) defined by (8.18) are orthogonal on the same set of spectral points but the measure 
acquires an additional linear factor 1 + x: 



^2 p s (l + x s )Q n (x s )Q m (x s ) = 0, if 



in 



Clearly, we have the case of a dense point spetrum on the interval [—1,1]. 

In the limit k = 1 we obtain the recurrence coefficients for the polynomials P n (x) 

1 (l-g n ) 2 (l-g 2n - 1 ) 2 (l + g n - 1 ) 2 

n 4 ( 1 + g 2n-l)2( 1 + g 2n)( 1 + g 2 n -2) 

, =n n-i/2 2( g +l)g"-(l-g)(l-g 2 ") 
" q (l + g 2 "- 1 )(l + g 2 "+ 1 ) 

where q = e~ 2w (of course, this q should not be confused with already introduced "elliptic" q by (3.11)). 

It is easily verified that these recurrence coefficients correspond to the Askey- Wilson orthogonal polynomials 

[4], [16] p n {z; a, b, c, d\q) with the parameters a = c = q 1 / 2 ,b = — d = 1. When max(\a\, \b\, \c\, \d\) < 1 then the 

Askey- Wilson polynomials are orthogonal on the interval [—1,1] with respect to a positive continuous measure 

[4], [16] 



/ 



1 w(x;a,b,c,d\q) , „ u , x 

p n {z; a, b, c, d\q)p m {z; a, b, c, d\q)dx = h n S n 

l v 1 — x z 
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where 

w(x-a b c dig) = M*;i)M*;-i)Ms;g 1/2 )M*;-g 1/2 ) 

' ' ' ' h(x; a)h(x;b)h(x; c)h(x; d) 

and 

oo 

h(x; a) = Y[ [1 - 2aa;<f + a V s ] 
Thus in our case the weight function becomes 

,_ fe(^;-g 1/2 ) _fr i + 2^ +1 / 2 + g 2s+1 ..... 

fc(a:;gV2) ~ 1 = 1 1 - 2xq*+V 2 + [ ' 

Expression (8.23) can be compared with the expression of the Jacobi dn-elliptic function in terms of the infinite 
product [24] 

We see that the weight function w(x) for the corresponding Askey- Wilson polynomials is proportional to the 
elliptic dn(u) function under the identification q — > q 2 , x — > cos(20). 

For the case of Q-polynomials defined by recurrence coefficients (8.16) we have in the limit k = 1 again the 
Askcy- Wilson polynomials p n (x; a, b, c, d\q) with the parameters a = c = q x / 2 ,b = l,d = —q. For the weight 
function we obtain 

which differs from the weight function (8.23) only by a linear factor 2(1 + x). 

It is interesting to note that despite the hyperbolic limit k — > 1 some "smell" of the elliptic functions remains 
in the expression for the weight function of the Askey- Wilson polynomials. Askey and Wilson themselves 
discussed intriguing cases when the weight function w(x) is expressible in terms of the Jacobi elliptic functions 
sn(it), cn(tt), dn(w) [4]. We see that the case of dn- weight function of the Askey- Wilson polynomials has a natural 
"elliptic" generalization leading to orthogonal polynomials with dense measure on the same interval [—1,1]. It 
would be interesting to search other possible elliptic generalizations of the Askey- Wilson polynomials. 



9. A general scheme leading to a dense point spetrum on the unit 
circle 

The examples of polynomials on the unit circle constructed in the previous sections can be generalized. In this 
section we propose a wide class of POC having a positive dense point spetrum on the unit circle. 
Let f(x) be a continuous function of a real variable x with three main properties: 

(i) f{x) is an even function f{—x) = fix) normalized by the condition /(0) = 1; 

(ii) function f(x) is periodic with some real period T: f{x + T) = fix); 
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(iii) the function f(x) has the Fourier expansion 

oc- 

/(*)= £ A n e 2 — / T , A_ n = A„, (9.1) 

n— — oo 

such that all the Fourier coefficients are nonnegative: A n > 0, n = 0, 1, 2, Moreover, we assume that there 

are infinity many positive coefficients A n > (otherwise the problem is trivial). Of course, we assume that the 
Fourier series converges everywhere on the real axis. 
We define the moments by the formula 

c n = f(wn), n = 0,l,2,..., (9.2) 

where w is an arbitrary positive parameter such that 

w ^ TM/N (9.3) 

for some integers M, N. Then it is obvious that c_„ = c„ and the moments c„, n = 0, 1, 2, . . . are all distinct: 
c n 7^ c m if 7i 7^ m. Using these moments c„ we can construct the Toeplitz determinants A„ and corresponding 
OPC by (1.2). 

Theorem 3 Under conditions (i)-(iii) the Toeplitz determinants A„ are all positive A„ > 0, or equivalently 
all reflection parameters satisfy the restriction — 1 < a n < 1. The polynomials P n {z) are orthogonal on the unit 
circle with respect to a positive dense point spetrum 

oo 

^ p s $„(z s )$ m (l/z s ) = /i„i„ (9.4) 

S — — oo 

where the orthogonality grid on the unit circle is 

z s = e 2mws ' T (9.5) 
and the discrete weights coincide with the Fourier coefficients 

Ps=A s , s = ±l,±2,... (9.6) 
Proof. From the Fourier expansion (9.1) we have the representation for the moments in the form 

OO OO r, 

c n = f(wn) = J2 A * z " = E P' z > = J z n Mz), n = 0, ±1, ±2, . . . 

S— — oo s— — oo ^ 

where z s is given by (9.5) and dfi(z) is the measure on the unit circle corresponding to the positive discrete 
weights p s > 0. Whence the polynomials P n {z) possess orthogonality property (9.4). The measure fj,(z) is 
positive and dense on the unit circle. Indeed, positivity follows from the property A n > for infinitely many 
values of n. On the other hand we see that in any small vicinity of the point z s on the unit circle there exists at 
least one point z s > from the set of spectral points z s due to condition (9.3). Thus the points z s form a Cantor 
dense set on the unit circle. This means that the corresponding point spectrum is dense. 
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Now from general properties of the polynomials orthogonal on the unit circle we conclude that all A„ are 
positive A n > and equivalcntly, that the reflection parameters satisfy the restriction — 1 < a n < 1. 

We thus see that starting from arbitrary even periodic continuous function f(x) having infinitely many 
positive Fourier coefficients A n we can construct corresponding OPC with a positive dense point spetrum on 
the unit circle. 

The Tocplitz determinants have the expression 



f 

/M 



i 



/W«-l)) 
/W«-2)) 



(9.7) 



/W«-l)) /W«-2)) •■• 1 

All the reflection coefficients of these polynomials will satisfy the restriction — 1 < a n < 1. Of course, if the 
function f(x) and its Fourier coefficients A n arc known explicitly, we have explicit orthogonality relation (9.4). 
However, explicit expression for the reflection coefficients a n can be obtained only for some exceptional cases. 

As far as we know the first (and so far the only) explicit example of such type was proposed by Magnus [17], 
[18] who considered the function f(x) defined as 

f{x) = {-\)W + \2{x-[x\) + \), 

where [x] means integer part, i.e. the greatest integer less than or equal to x. Equivalently, the function f(x) 
can be defined as f(x) = 1 — 2\x\ for —1 < x < 1 and then continued for the whole real axis by 2-periodicity. 

Clearly the function f(x) is even /(— x) — f(x) and periodic with period 2: f(x + 2) = f(x). Moreover, 
there is an obvious property f(x + 1) = —f(x). If two numbers x and y have the same integer part [a;] = [y] 
then obviously 



m - f(y) 



x-y 



= -(-l) [ 



(9.8) 



The function f(x) has the Fourier expansion 



g7rza;(2s+l) 

(2s +1) 2 



Thus A n = 4/(7r 2 n 2 ) if n is odd and A n = if n is even. We see that all A n are nonnegative, hence conditions 
for existence of a positive dense point spetrum are fulfilled and corresponding polynomials are orthogonal on 
the unit circle 

oo 

(2s + l)- 2 <P n (z s )<S> m (l/z s ) = h n 5 nm (9.9) 

S— — OQ 

with some positive normalization coefficients h n . The spectral points on the unit circle are 



i7r(2s+l)w 
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Magnus showed [18] that the polynomials § n {z) (and corresponding reflected parameters a n ) can be found 
explicitly. This explicit expression is closely related with the expansion of the irrational number w to the 
continued fraction 

w = q + , (9.10) 

1 

9i + 



<72 + 



q 3 + ... 

Introduce corresponding convergents Pi/Qi, i = 0,1,2,..., where P = [u>],Qo = 1 is "zero-step" rational 
(integer in this case) approximation of the number w and the rational number P n / Q n is obtained when truncating 
the continued fraction (9.10) at the n-th step: 

Pn 1 

q- = qo + 



1 

q 2 + ... + — 

Qn 



The n-th convergent P n /Q n provides the best rational approximation of the irrational number in the following 
sense: 

\Q n w - P n \ < \yw - x\ 

for all possible integers x,y with the restriction y < Q n . 

For any given positive integer n we define a sequence of pairs of nonnegativc integers (n^, rrii) in the following 
way. Let the pair (rij, mi) provides the z-th order best rational approximation of the number w under condition 
1 < rii < n. In more details this means the following. Define positive (irrational) numbers yi = \wrii —mi\. The 
pairs (rij, mj), i = 1, 2, . . . , n correspond to the ordering < y\ < y 2 < y 3 < ■ ■ ■ < y n , where 1 < rij < n, i.e. the 
pair (ni,rni) yields the best rational approximation to the number w (with denominators not exceeding n), the 
pair (712,7712) yields then second (next) best approximation etc. It is well known that the best approximation 
(ni,TOi) is given by an appropriate continued fraction of the number w. mi = ri\ = Qk, where Qk < n. 
The second approximation (n 2 ,m 2 ) as well all other approximations (n^rrij) can be easily expressed in terms 
of so-called intermediate fractions (see e.g. [15]). Now we can always present the polynomial $„(z) in the 
following form 

$„(*) = z n + Gipz n - ni + G^z n ~ n2 + ... G^z 11 - 11 * + . . . + Gl n) (9.11) 

(s) 

with some coefficients G n , where m, n 2 , . . . correspond to the succeeding best approximations of the number 
w. (Clearly, the set m, n 2 1 ■ ■ ■ n n is a permutation of the set 1, 2, . . . , n, i.e. n, 7^ nk if i 7^ k due to irrationality 
of w). It can be easily showed (see [18] for details) that only 3 first terms survive in (9.11), i.e. we have 

= z n + G^z n - ni + G^z n - n2 (9.12) 

where the coefficients Gn \ Gn" 1 can be found explicitly if {n\,m{) and (n 2 , m 2 ) are known [18]. 
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Note that in general we have a n = for "almost all" n. The nonzero coefficients a n appear only if for some 
n we have m = n or ni = n. We thus need information about the best approximations of the number w. 
Unfortunately for almost all irrational numbers w it is impossible to present some "explicit" formula for n\,ni 
(expressing ni,n 2 as an analytic function of n). 

In contrast, in our examples of cn- and dn- elliptic polynomials we have purely explicit presentation for 
polynomials & n (z), moments c„ and reflection parameters a n . 

Is it possible to construct other explicit examples of the polynomials $ n ( z ) corresponding to functions f(x) 
satisfying properties (i)-(iii)? 



10. Remarks concerning "elliptic hyper geometric functions" 

Consider cn or dn polynomials § n {z) on the unit circle as some "elliptic hypergeometric functions" depending 
on the argument z and discrete parameter n. The parameter w plays the role of the "deformation" parameter: 
when w — Q we have from (3.1) 

$„(*) = (*-l) n (10.1) 
for both cn and dn polynomials. Thus in the limit «) = 0we have the simplest hypergeometric function 

lim $„(z) = (-1)" iF (-n; z) = (z - 1)". 

w— >0 

This allows us to consider the functions ^^\z) and &lf\z) as "elliptic deformations" of the hypergeometric 
function iF (— n;z). 

Note that from the explicit expression (3.1) it is sen that functions <&^n\z) and $^(z) do not belong to a 
class of elliptic hypergeometric functions 12V11 introduced by Frenkel and Turaev [10] and studied intensively 
during last years (see e.g. [21], [11] and many others). In particular, these functions perhaps do not possess 
simple transformation formulas like Bailey transform [21]. Nevertheless, functions <^^\z) and <&\®\z) have 
several nice properties which are very close to the classical ones. In particular the recurrence relation (1.12) can 
be considered as a contiguous relation for 3 elliptic hypergeometric functions $„(z) with parameters n, nil. The 
transformation properties (5.5) with respect to the operator £ can be considered as an elliptic analogue of the 
derivatives of the hypergeometric functions. Finally, the functions $l C ^(z) and ^^\z) possess nice summation 
formulas for two values of the argument z = ±1. Indeed, from (1.21) we have, e.g. for the cn-polynomials 

$ 2n } (!) = - cn(w))(l + dn(2w)) ... (1 - cn(w(2n - 1))(1 + dn(2wn) 

and 

$2^+1(1) = C 1 - cn(u?))(l + An{2w)) . . . (1 + dn(2um)(l - cn(w(2n + 1)) 

(similar expressions hold for $^(±1)). 

Thus there exists a much simpler class of " elliptic hypergeometric functions" with nice properties similar to 
the classical ones. 
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